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1. Introduction 

Determining the spectrum of BPS states in four-dimensional gauge theories and string vacua 
with Jf = 2 supersymmetries has been an active topic of research since the seminal work of Seiberg 
and Witten [1]. Indeed, it affords a window into strongly coupled phenomena triggered by non- 
perturbative states, such as monopole confinement and condensation in softly broken N = 2 theories 
or conifold transitions between apparently disconnected Calabi-Yau (CY) vacua, and allows for 
precision tests of non-perturbative dualities. Notably, it remains a challenge to determine the exact 
degeneracies of BPS black holes in jV = 2 string vacua, and study deviations from the classical 
Bekenstein-Hawking entropy at finite electromagnetic charges, to the same level of precision as 
has been achieved for jV = 4 or jY = 8 string vacua (see e.g. [2, 3, 4, 5, 6] for the current state of 
the art). 

While the BPS spectrum is largely insensitive to variations of parameters or coupling con- 
stants, as BPS states lie in short multiplets and can only desaturate the BPS bound by combining 
with another BPS state of opposite parity, the BPS index £l{y,t) = Tr(— l) 2j3 is in general not 
constant over the space SS of parameters t : the reason is that along certain real-codimension one 
walls W in SB, known as walls of marginal stability, the point-particle spectrum meets with the 
continuum of multi-particle states, signaling that a BPS state with charge y is about to decay into 
more elementary BPS constituents with charge a, with the same total charge y = £,a, [l] 1 . Across 
the wall, the total BPS index £l{y;t) jumps by a quantity which is equal to the BPS index of the 
bound state, which contributes on one side of W but not on the other. While this wall-crossing phe- 
nomenon complicates the analysis of the BPS spectrum, it turns out to be governed by a completely 
universal formula, which holds both in gauge theories and string vacua. This formula first appeared 
in various guises in the mathematics literature on Donaldson- Thomas (DT) invariants [8, 9], and 
has received various physics proofs since then [2, 10, 11, 12, 13, 14, 15, 16]. 

In the first lecture, we shall outline a physically transparent derivation of this wall-crossing 
formula [15], which hinges on a general solution of the supersymmetric quantum mechanics of 
multi-centered solitons. While being equivalent to the original wall-crossing formulae of [8, 9], its 
structure is completely different and has already led to a number of new mathematical insights on 
DT invariants [17, 18] (see [19] for an earlier review of wall-crossing, with a different emphasis). 

One implication of wall-crossing is that at least some of the BPS states are not elementary, 
being liable to decay into more elementary constituents. There are reasons to believe, however, 
that there might exist a subset of BPS states which are absolutely stable throughout the parameter 
space £$. This is certainly the case in the simplest jV = 2 gauge theory considered in [1, 20], 
where the monopole and dyon are the only BPS states in the strong coupling chamber, while all 
BPS states in the weak coupling chamber arise as bound states of these two (see Figure 1). It is 
also clear in string theory vacua that BPS states represented by single-centered black holes cannot 
decay (except perhaps at the boundary between two basins of attraction [21]). Thus, it should be 
possible to reconstruct the total BPS index £l{y\t) from the BPS degeneracies £2s(y) associated to 
these 'single-centered' or 'elementary' constituents, by solving the quantum mechanics of a BPS 
bound state of n single-centered constituents with charges {a,}. Since the charges are unrestricted, 
the quantum mechanics is in general more involved that in the case relevant for wall-crossing, 



'A similar phenomenon was first encountered for two-dimensional theories with (2,2) supersymmetry in [7]. 



Wall-crossing, multi-centered black holes, and quiver invariants Boris Pioline 



(2,0) 

(2n,l) 

(2n+2,-l) 




Figure 1: BPS spectrum in jV = 2,D = 4 SYM theory with SU(2) gauge group and no flavor. In the strong 
coupling chamber, the only stable BPS states in the strong coupling chamber are the monopole and dyon, 
with charge (0,-1) and (2,-1). All states in the weak coupling chamber outside the wall arise as bound 
states of the monopole and dyon. 

as the phase space of multi-centered black holes can have non-compact regions where some of 
the centers can approach each other to arbitrarily small distance. In Lecture 2 I shall present a 
general prescription for computing the index of such bound states, which is based on the physical 
hypothesis that these 'scaling solutions', which classically have zero angular momentum, carry the 
smallest possible angular momentum at the quantum level [22]. 

While it is at present impractical to compute the total BPS index £l(y, t) or single-centered BPS 
index Qs(y) in general Calabi-Yau string vacua, as it would require determining the complete set 
of DT invariants of the CY threefold 36 ', it is sometimes possible to realize a class of BPS states 
as bound states of certain D-branes wrapped on supersymmetric cycles in 3£ , and described by 
some matrix quantum mechanics described by a certain quiver [23, 24]. Mathematically, the BPS 
bound states are represented by cohomology classes on the moduli space of semi-stable quiver 
representations. This provide a concrete framework to test the afore-mentioned hypothesis, and 
identify the elementary constituents as special cohomology classes. A recent study of three-node 
quivers with a loop [25] and similar necklace quivers [26] suggest that these elementary states can 
be understood as middle cohomology states, or 'pure Higgs' states, which turn out to be remarkably 
robust under wall-crossing. In Lecture 3 we shall review the techniques used to analyze these 
quivers and use our general prescription to identify the 'pure-Higgs' states in a more precise way 
[27]. 

2. Wall-crossing from multi-center black hole quantum mechanics 

In this first lecture, we discuss wall-crossing phenomena in four-dimensional string vacua 
described by JV = 2 supergravity, or in jV = 2 gauge theories described by an Abelian gauge 
theory on the Coulomb branch. In both cases, the jump of the index across the wall is universally 
determined by the quantum mechanics of multi-centered solitons, which is identical for black holes 
and dyons. 

2.1 Generalities 

Let r = T e ©r m be the lattice of electric and magnetic charges (q\,p A ), with symplectic 
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(Dirac-Schwinger- Zwanziger) integer pairing 

(r,/)=w' A -WAeZ (2-1) 

BPS states preserve 4 out of 8 supercharges, and saturate the bound M > \Zy(t)\ where the central 
charge Zy{t) = (Y(t),y) is a linear functional of the electromagnetic charges [28]. The index 



1 

2 



a{y;t) = --Tr^ w (-l) 2 ^(27 3 ) 2 = Tr W-l)^ (2.2) 



(where My(t) is the Hilbert space of one-particle states with charge y € T in the vacuum with vector 
moduli t £ 38, and the prime denotes the removal of the center of motion degrees of freedom) 
receives contributions from short multiplets only. In JV = 2 gauge theories (but not in jV = 2 
string vacua), there is an additional SU(2)r symmetry, which allows to defined the protected spin 
character (PSC) CL(Y,t,y) = Tr(-l) 2j3 y 2( - h+J ^ [13]. The 'non exotics' conjecture [13] asserts that 
all states have Ij = 0, and the PSC coincides with the refined index, defined by 

£l( r ,t,y) = Tr^, {t) (-l) 2J iy 2J i. (2.3) 

In ,jY = 2 SUGRA/string vacua, one can still define Q.(y;t,y) by (2.3) but it is no longer protected, 
hence it could get contributions from non-BPS states and depend on HM moduli. Nevertheless, it 
will be a very useful device to compute the protected index (2.2), obtained by setting y = 1 in (2.3). 
The BPS index 0.(y,t) is locally constant, but it may jump on codimension-one loci 

W{YuY2) = {t € 38 : Zy.it) /Z^t) £ R+} (2.4) 

in vector-multiplet moduli space, called 'walls of marginal stability'. The basic mechanism is that 
some of the BPS states with charge y = My\ +Nj2 are bound states of more elementary BPS states 
with charge a,- = M,yi +NiJz with Y,j{Mi,Ni) = (M,N), which exist only on one side of the wall 
W(Yi,Y2.). Several remarks are in order: 

i) for a given wall, we can always choose the basis vectors y , j2 to be primitive (i.e. such that 
Yi/d $. r for all integers d > 1), and such that Q.(My\ +NYi) has support on the positive 
cone r = Ny +Ny2 (corresponding to particles) and its opposite — T (corresponding to anti- 
particules) [14]. The constituents lie in f , so only a finite number of bound states can occur. 

ii) The index of states with y ^ Zy +Zy? is constant across the wall W{Y\,Yi), and so are£2(y) 
and Q(ya). 

iii) The index Q.{My\ +NY2) may contain contributions from bound states of constituents with 
charges lying outside the sublattice Zy + Zy2 C T, but those are insensitive to the wall. 

iv) The states {a,} do not form a BPS bound state on the wrong side of the wall, but one cannot 
exclude the possibility that they could form a non-BPS bound state. Such states do not 
contribute to AH(y), nor to the PSC, but in string vacua they could contribute to the refined 
index A£l(y,y). We shall work under the assumption that this possibility does not occur. 
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2.2 Primitive wall-crossing 

As a warm-up we consider the fate of BPS states with charge 7 = j\ + 72 across the wall 
W{j\ , 72). For (71 , 72) 7^ 0, there exists a two-centered BPS solution of charge 7 = 71 + 72, angular 
momentum J = \ (71 , 72) r/|r| of the low energy effective action, with separation given by [29] 

m = (71,72)1^(71+72)1 (2 _ 

1 ' 23[Z(7i)Z( 72 )] l ■ 

The solution exists only on one side of the wall. As t approaches the wall, the distance \r\ diverges 
and the bound state decays into its constituents j\ and 72. Near the wall, the two centers can 
be treated as pointlike charged particles with Cl(yi) internal degrees of freedom, interacting via 
exchange of massless particles, including the usual Coulomb and Lorentz forces. The classical 
BPS phase space ^#2 for the two-particle system is the two-sphere, with symplectic form [24] 

(0 = \yn sin dddd(j) , 7i 2 = (7i,7>) (2-6) 

such that rotations d^ are generated by the Hamiltonian 2 J$ = \j\i cos 6. The geometric quantiza- 
tion of ^#2 leads to I712I states transforming as a spin j = j(\Yn\ — 1) multiplet under rotations. 
Near the wall, these conflgurational degrees of freedom decouple from the internal degrees of free- 
dom of the two particles. The index of the two-particle bound state is then ( — l) 7l2+1 7i2 x ^(/i) x 
£2(72). These are the only bound states of charge 7 = 71 + 72 which (dis)appear across the wall, 
thus the jump of the index is given by the primitive wall-crossing formula [2] 

An(7 1 + ^) = (-l)^ +1 712^(70^(72) . (2.7) 

Similarly, the variation of the refined index is given by [10, 30] 

An(7 1 + 7 2 ;y) = l y _\ jy — ft(>i;y)fl()4;y) , (2.8) 

which reduces to the previous formula at y = 1 . 

2.3 The quantum mechanics of multi-centered black holes 

We shall now generalize this reasoning to compute AQ.(Mji +./V72) for general (M,N). For this 
we need a description of the BPS phase space of multi-centered black holes (or dyons). The most 

general stationary, n-centered BPS solution of JV = 2 SUGRA with total charge y=a\-\ \-a„ 

is [29, 31] 

ds 2 = -e 2U {At + £f) 2 + e~ 2U dr 2 (2.9) 

where the scale factor U (r) and the scalar fields t (?) are determined by solving 

2e- u ^3 [e-'^Y (t(r))] = J3 + £ -^— . (2.10) 



Recall that the Hamiltonian (or moment map) of a vector field K preserving the symplectic form ft) is the function 
fl such that l K a> = d/x. The additive ambiguity is uniquely fixed by the condition {[Ik,Hk'} = IMk k'] m me case °f a 
simple Lie group action. 
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The one-form &/ and the electromagnetic fields are uniquely determined, provided the locations r, 
satisfy the conditions 



Y ( ffl^A = Ci . (2.11) 






In these equations, we have denoted 

0=argZ 7 , , J8 = 2Im[«T i *y (*„„)] , a = 2Im(e^Z«,) . (2.12) 

Eq. (2.11) enforces n — 1 conditions on 3rc locations r,-. Modding out by translations in IR 3 , the 
space of solutions is a 2rc — 2 dimensional space which we denote by ^„({a,},{c,}). For the 
solution to be regular, one must also ensure 



2C ^) = -S /J + V-^L- >0, Vrei 3 , (2.13) 




e 

where 5(j8) is the Bekenstein-Hawking entropy functional [31]. While this condition may in gen- 
eral remove some connected components in ^ n , it appears to be automatically satisfied when all 
charges a r - lie in a two-dimensional lattice Mj\ + Nyi with MN > 0, the case relevant for wall- 
crossing. In this case, the constants c, are uniquely determined by c,- = A£^,-(CU;,GJ/}, up to an 
overall scale A which vanishes on the wall. For =yf = 2 gauge theories, the most general multi- 
centered solution of the low energy Abelian gauge theory on the Coulomb branch has not been 
worked out, but it is widely believed that solutions are similar to the one above, with the same 
equilibrium condition (2.11) on the center (but without the constraint (2.13)) [32]. 

When the centers are well-separated, the dynamics of the bound state is described by JV = 4 
quantum mechanics with 3« bosonic coordinates, An fermionic coordinates with Lagrangian [24, 
33, 34, 32] 

& = YjWift)] 2 + L^A + L MO' + • • • (2- 14) 

i i ij 

W t = £ ±JLJL- Ci , A . = £(«., a 7 .)A Dirac (r r7 ) , (2.15) 

where r,y = r,- — rj, and Ach-ac^) is the gauge potential of a unit charge Dirac magnetic monopole 
sitting at r = 0. Factoring out the center of mass motion, the classical ground state dynamics 
is then first order quantum mechanics on the BPS phase space ./#„({a;},{c,}) = {W, = 0}/R , 
equipped with a symplectic form ft) uniquely determined by the requirement that the moment map 
of the action of SO(3) rotations on j$ n coincides with the total angular momentum / earned by the 
configuration [24, 33] 

ft) = jJ/Cka^dADfcacOT/) i ^=jE(<%.0;)- • ( 2 - 16 ) 

The geometric quantization of .M n produces a BPS Hilbert space given by zero-modes of the natural 
Dirac operator on (^#„,co). The refined index (or PSC in gauge theory setting) is equal to the 
equivariant Dirac index [22, 34, 32] 

gCotiomb({0Ci},{ci};y) = Tr Ke rZ) + (-v) 273 - Tr Ke rD_ {-yf h ■ (2.17) 
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For reasons which will become clear in Lecture 3, we refer to (2.17) as the Coulomb index of the 
n-centered configuration. In the limit where all DSZ products a,-/ = (a,-, Of/) are scaled to infinity, 
the Coulomb index reduces to the equivariant symplectic volume 

gctesdOiUcihy) = \' 1( 7vr / co'^y 2 ' 3 ■ (2.18) 

As we shall see, the index can be computed exactly for any number of centers using localization 
with respect to rotations J3 around a fixed axis. 

2.4 Localization and the Coulomb index 

Away from walls of marginal stability, defined by hyperplanes £, GA c\ = in the space of the 
parameters c* (where A is any proper subset of {1, . . . ,n}), the relative distances |r, ; | are bounded 
both from above and from below, < r < |r,/| < R so that the BPS phase space j% n is compact. 
Moreover, the action of ^3 has only isolated fixed points, namely collinear configurations along the 
z axis, subject to the one-dimensional reduction of (2.1 1), 

n (a- a-) 1 

£ H^r =c {) h = - £(«,-, a } ) sign(z ; - - zi) . (2.19) 

These are n— 1 conditions on n — 1 locations Zi~Zi, hence typically admit isolated solutions. Thus, 
both gciass and gcouiomb can be evaluated by localization with respect to the rotations ^3 along the z, 
axis, using the Duistermatt-Heckman and Atiyah-Bott-Lefschetz fixed point formulae, respectively 
[15, 22]. For (2.17) this gives 

where L is the matrix of the action of 73 on the holomorphic tangent space around the fixed point. 
With some work, the determinant around a fixed point (2.19) evaluates to (y — l/y) n ~ [ times a sign 
s(p) = — sign(detH /// ) where W" is the Hessian of the superpotential 

W({zt}) = -^{ai,aj} sign[z j-Zi] \a\zj-zt\ -J^c/z,- , (2.21) 

><j > 

whose critical points reproduce the constraints (2.19). After the dust settles, one finds that the 
refined Coulomb index is given by [15, 22] 

gCoulomb({^},{c«};j) = (-1) L< ^" l £, ( ) ^, <J ( a ,,« j )signfa- z ,) (2 22) 

(y-y 7 p 

where the sum runs over all solutions of (2.19). The contribution of each fixed point is singular at 
y = 1, but the sum over fixed points is guaranteed to produce a symmetric polynomial in y and l/y, 
as long as j% n is compact. In the simplest two-centered case with c\ (0:1,0:2) > 0, the fixed points 
on the BPS phase ^2 = S 2 , lie at the north and south pole, with 2/3 = ±(0fi , Ola), hence 

v («l,«2> _ V -<«1, Oi) 

Scouiambttai.ofc}) =(-l) a - +lJ f- = Tr , fi . (2.23) 

y—l/y j=2({ui-ai)-\) 
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If the sign of c\ is opposite, the BPS phase space is empty and gcouiomt>({«i , 0^2}) vanishes. Setting 
CCi = 71 , 0J2 = 72 and multiplying by the degeneracies Q.(ji;y) carried by each center, one recovers 
(2.8). 

For more than 2 centers, the enumeration of solutions to (2.19) becomes rapidly impractical, 
even on a computer. It is worth however noting that the contribution of a given ordering of the cen- 
ters to (2.22) does not depend on the details of the solution, only on the Morse index of the critical 
point of W, and so is invariant under deformations of the DSZ products CCtj and constants a, away 
from certain singular loci. This robustness was exploited recently to give an inductive algorithm for 
computing the Coulomb index, which is much more efficient than looking for numerical solutions 
to (2.19) [35]. In the case relevant for wall-crossing, this algorithm gives a fully explicit result, 

8cou l orn b ({au---a n y,{c u ---c n y,y)=(-\) n - 1+ ^ a Hy-y- 1 )- n+i 

"- 1 / * \ _j (2 24) 

LEI «*(*),*(*+!) L C *W (-l) E *=i (- a -«.^))y E ^ a -('>w , 

<T k=\ \ i=\ J 

where &(x) is the Heaviside function, equal to 1 if x > and otherwise, a runs over all permuta- 
tions of { 1 , . . . , n } , and we denote a i; - = ( a, , ccj ) . 

2.5 Statistics and the Coulomb branch wall-crossing formula 

Having computed the index of the quantum mechanics of n centers, we can apply the same 
logic as before and tentatively write the jump in the total index as a sum of the indices of all bound 
states that decay across the wall, 

Aa( 7; j) = £ £ gdaiYrflYLi^fay) < 2 - 25) 

n>2 {a,,...a„}ef 
7=a,H — ha„ 

where 7 = Mj\ + ./V72, a r - = Mij\ + Niji, and Q. + ((Xf,y) is the refined index on the side where the 
bound state does not exist. This is almost right, but it overlooks the issue of statistics. If the centers 
were classical, distinguishable objects, Maxwell-Boltzmann statistics would require to weight each 
(unordered) contribution {a,} with a symmetry factor l/|Aut({a,})|, where |Aut({a,-})| is the 
order of the subgroup of the permutation group leaving {a,} invariant. Instead, the centers are 
quantum objects, with Bose statistics if 0(a,) > 0, or Fermi statistics if Q(a,) < 0. One can show 
that the Maxwell-Boltzmann prescription nevertheless works, provided one replaces everywhere 
the index £l(y) with its rational counterpart [15] 

A(roO - 1 J7-vV (yA/y) ' ^ (7) = E T^ {r/d) ■ (2 " 26) 

d\y \J J ) d\y 

To see why this can work, consider for example AQ.(ji + 272): it receives contributions from bound 
states {71 +72,72}, {71,272}, {71,72,72}- The last, 3-centered case can be treated as a two-body 
problem, since the two centers with charge 72 have zero DSZ product hence do not interact. Taking 
into account Bose-Fermi statistics for {71 , 72, 72}, one finds, setting y = 1 for simplicity, 

AQ(7 1 +272)=(-lp7i2^ + (72)^ + (7i + 72)+27i2n + (272)a + (7 1 ) 
+ 57i 2 n + (72)(7 12 ^ + (72) + l)^ + (7i) • 
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The last contribution violates the naive charge conservation implied by (2.25), but is forced by 
Bose-Fermi statistics. Rewriting this formula in terms of the rational invariant A(27 2 ) = £2(2/2) + 
^Q.(Yi), one finds that naive charge conservation is restored, with the expected | Boltzmann sym- 
metry factor, 

AA(7 1 +27 2 )=(-l)^7 12 A+(7 2 )A + (7 1 + 7 2 )+27 12 A + (27 2 )A + ( 7l ) 

+ 2 -[yi 2 6 + (72)] 2 6 + (yi)- 



(2.28) 



A general justification for the replacement Q. — > Q. can be found in [15]. 

After these amendments, we arrive at the Coulomb branch wall-crossing formula, 

«>2 {a 1 ,...a„}er ' VL l > n 

Remarkably, this formula can be shown to agree with the wall-crossing formulae established in 
the mathematical literature on Donaldson- Thomas invariants for the derived category of coherent 
sheaves on a Calabi-Yau threefold. The agreement is far from trivial and can be established by 
induction on the number of centers [36]. As a simple example, the formula (2.29) with (71 , 72) = 2 
reproduces the weak coupling spectrum of pure SU(2) gauge theory, from the knowledge of the 
strong coupling spectrum (see Figure 1) [11]. 

3. The total index from single-centered degeneracies 

The wall-crossing phenomenon discussed in Lecture 1 shows that the BPS states contributing 
to the total refined index (or PSC) Q.(y;t;y) are in general not elementary, but can be decomposed 
into more elementary constituents. This suggests the possibility that there could exist a set of 
elementary, absolutely stable BPS states, with moduli-independent indices Q.s(y;y), such that any 
other BPS state would arise as a bound state of those. This is indeed realized in pure SU(2) gauge 
theory, where all states arise as bound states of the monopole and dyon. More generally, it happens 
in any jV = 2 gauge theory which has a strong coupling chamber with a finite number of BPS 
states. In JV = 2 string vacua, it is clear that BPS states represented by single centered black 
holes cannot decay and so should be considered as elementary. In this lecture, I shall address the 
question of computing the total index Q.{y;t;y) (not only its variation) in terms of indices Q.^{y;y) 
associated to elementary constituents, without committing to the nature of these constituents. 

According to the physical picture explained in Lecture 1 , it is natural to expect that the total 
index Q.(y;t;y), at any point in moduli space, should be a sum of contributions of all possible bound 
states of elementary constituents with charge a, such that £ a, = 7 Assuming that the BPS Hilbert 
space associated to each bound state factorizes into configurational degrees of freedom, described 
by the geometric quantization of the BPS phase space ^#„({a,},{c,}), and internal degrees of 
freedom Q.s(cCj;y) associated to each center, and taking into account Bose-Fermi statistics, one 
may tentatively write 

y=aiH ya„ 
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where the rational invariants £ls((Xi,y) are defined in terms of the integer invariants D.$(cCi',y) as i n 
(2.26). While similar in spirit with the Coulomb branch wall-crossing formula (2.29), this formula 
differs in several important ways: 

i) it includes a contribution from single centered black holes, with n = 1 ; 

ii) the charges a, run over all possible vectors in F with &s{cCi) / 0, rather than a two-dimensional 
sublattice f ; 

iii) the Coulomb index gcouiomb({cd},{ci},y) depends on the moduli at infinity t through the 
constants c,-; in the context of ,jV = 2 supergravity it should only count states consistent with 
the regularity condition (2.13). 

The fact that the charges of the constituents are no longer restricted to the positive cone in a two- 
dimensional sublattice has drastic consequences, as we now discuss. 

3.1 Scaling solutions 

Unlike the case relevant for wall-crossing, the BPS phase space ^#„({a,},{c,}) for generic 
choices of the DSZ product is in general non-compact, even away from walls of marginal stability. 
To see this, consider the three-centered case with DSZ products a\ 2 = a, a 2 j, = b, (X31 = c satisfying 
triangular inequalities < a < b + c, etc. It is easy to check that the equations (2. 11), 

a c b a c b 

=ci, = c 2 , =c 3 , (3.2) 

m ri3 r 2 3 rn m r 23 

admit solutions where the three centers can approach each other arbitrarily, r\ 2 ~ as, r 23 ~ &£, H3 ~ 
c£ with £ — > [2, 37]. These configurations, known as 'scaling solutions', exist for any choice of 
c\,c 2 ,C3 such that c\ +c 2 + C3 = 0, and carry arbitrarily small angular momentum J 1 ~ e 2 . For 
c\,c 2 > 0, a naive application of the Atiyah-Bott fixed point formula (2.22) gives 

( _\\a +b+c (y a+h ~ c + y- a ~ h + c ") 
^Coulomb — 1 -, 1 no > w--V 

{y-^/yr 

where the two terms arise from the regular collinear configurations (123) and (321). This is not a 
polynomial in y, in particular it is singular as y — > 1 . An explicit computation shows that the equiv- 
ariant volume of ^#3 is nevertheless finite, and agrees with the Duistermat-Heckman formula after 
including by hand an additional fixed point contribution from the scaling boundary of ^#3, with 
effective angular momentum Jj = [22]. Similarly, (3.3) can be turned into a Laurent polynomial 
for any a, b, c by adding by hand a fixed point contribution with Jt, = (when a + b + c is even) or 
J 3 = i (when a + b + c is odd): 



i\a+b+c I yO+b— c 1 y— a— b+c . 



<?Coulomb 



2 , a + b + c even 
y + 1 jy , a + b + c odd 



(3.4) 



{y-l/yY 

More generally, scaling regions in ^#„({a,}, {c,}) arise whenever there exist a subset A and vectors 
fi eM 3 ,/EAsuchthat 

Vie A, Y ^ = . (3.5) 

jeA\ r 'J\ 
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Such solutions exist, when they do, in all chambers in the space of the parameters c,'s, and have 
vanishing angular momentum J 1 as the centers approach each other. Despite these non-compact 
regions, one can argue that the symplectic volume is finite, and it is plausible that ^#„ admits a 
compactification such that the Atiyah-Bott Lefschetz formula still applies (after allowing for non- 
isolated fixed points at the boundary). Indeed, it is possible to formulate a 'minimal modification 
hypothesis' which converts the Coulomb index gcouiomb computed from (2.22) into a symmetric 
Laurent polynomial gcouiomb > by adding by hand suitable fixed point contributions with the mini- 
mum possible angular momentum J^ [22]. Unfortunately, this is not good enough to ensure that 
the answer from (3. 1) will be a symmetric Laurent polynomial, as we now discuss. 

3.2 The Coulomb branch formula for the total index 

Even if the Coulomb indices gCouiomb({c^'},{ci},;y) were Laurent polynomials, the tentative 
formula (3. 1) would not necessarily produce a Laurent polynomial for Q.(y;t;y), due to the fact that 
the rational invariants Cls((Xi,y) are rational functions of y when a,- is not primitive. It is possible 
to give a minimal prescription that cures this problem by introducing new 'scaling' contributions 
as follows [22] 



n/„. \ V V gCoulomb({«<},{Q},>') 

<W)= |,eL |a-««»I 



n 



_y-\/y 
The new contribution £2 sca n ng is determined recursively in terms of Q.$ by 



(3.6) 



Ii I mi ^-%) [^(oa/mr^+a^^/m^)] 



n sca i ing (a;j) = £ H({p i };{m i };y)Ha s iP i ;y m ') (3.7) 

{fter}.{m ; ez} i 

where the rational function H({fij};{mi};y) is in turn determined uniquely by the 'minimal modi- 
fication hypothesis', 

i) H is symmetric under y — > \jy and vanishes as y — > 0, 

ii) the coefficient of floret' (A'./"') in the expression for £l(E,-wi, , j3,-;;y) obtained using (3.6) must 
be a Laurent polynomial in y. 

It can be checked that the condition ii) exactly fixes the free parameters allowed by i). As an 
example, for the 3-center configuration discussed in §3.1, 

, . | — 2(y-y~ 1 )~ 2 ,a + £> + ceven 

#({yi,72,73};{i,i,i};)0 = r * ' 2 , (3-8) 

[(y+y )(y-y ) ,a+b+codd 

so the prescription r(3.6) reduces to the replacement gcouiomb — > ^Couiomg. but this is not so in 
general. The Mathematica package developed in [35] provides a convenient implementation of this 
prescription. 

The Coulomb branch formula (3.6) conjecturally expresses the set of indices Q.(y;t,y) in terms 
of a new set of indices Q.^{y;y). One might ask, what have we gained ? First, unlike the total 
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index Q.{y;t;y), the single-centered invariants Q.v,(y;y) no longer depend on the moduli. This 
may be a valuable property in enforcing invariance under various dualities, since these dualities 
typically act both on the charge y and moduli t. Second, in j¥ = 2 string vacua we expect that 
Q.s(y;y) counts micro-states of single-centered BPS black holes, withAdS2 near-horizon geometry. 
Supersymmetry requires these BPS microstates to lie in singlets of the rotation group [38], thus 
their refined index Q.s(y,y) should be independent of v. If this hypothesis is correct, then the single- 
centered invariants Hs(y) provide a concise parametrization of the BPS spectrum in all chambers 
of moduli space, which can be directly compared with the path integral of quantum gravity in a 
single AdSj throat, also known as the quantum entropy function [39]. In the next lecture, we shall 
test this hypothesis in the case of quiver quantum mechanics, which provides a model for D-brane 
bound states. 

4. Applications to quivers 

In this last lecture, we shall apply the Coulomb branch formula (3.6) in the context of quiver 
quantum mechanics (or, in mathematical parlance, the cohomology of quiver moduli spaces), where 
the refined index Q.(y;t;y) is unambiguously defined and often computable. Quiver quantum me- 
chanics on the Higgs branch provides an accurate description of D-brane bound states near loci 
in moduli space where the central charges of the constituents nearly align [24]. On the Coulomb 
branch, it also reproduces the quantum mechanics of multi-centered solitons discussed in §2.3. 
This physical duality has very interesting mathematical consequences for the cohomology of quiver 
moduli spaces, as we shall explain. 

4.1 Quiver quantum mechanics 

Quiver quantum mechanics is a particular type of matrix quantum mechanics with jV = 4 
supercharges, whose matter content is encoded in a quiver, 3 , i.e. a set of nodes and arrows (see 
e.g. (4.6) or (4.15) below). Each node £ = 1 .. .K represents a U(Ne) vector multiplet if^D[), each 
arrow represents a chiral multiplet (p^j: in (N(,Nk) representation of U(Ne) x U(Nk)- The set {Afy} 
is called the dimension vector. In addition, one must specify Fayet-Iliopoulos terms a such that 
Y^Nice = 0, and (in presence of closed oriented loops) a gauge invariant superpotential W(<j>k,i). In 
applications to D-brane bound states, each node represents N( coinciding D-branes of charge ye and 
the number of arrows from node £ to k is given by the DSZ product 4 Ytk- The FI term a depends 
on vector multiplet (VM) moduli, while the coefficients of W in general depend on hypermultiplet 
(HM) moduli. 

On the Coulomb branch, the chiral multiplets and off-diagonal components of the vector mul- 
tiplets are massive and can be integrated out. The diagonal components of the vector multiplets 
re then represent the locations of N( Abelian D-branes carrying charge yi, subject to the D-term 
constraints 

W, £_% T = c,. (4.1) 

Quivers also describe SUSY gauge theories in higher dimension, but here we focus on D = 1 

In general, the DSZ product Jtk gives the net number of arrows from node I to k, counted with orientation, but we 

shall always assume that the superpotential W is generic such that chiral multiplets associated to pairs of arrows with 

opposite direction are massive and can be integrated out. 
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Remarkably, this reproduces the equations (2.11) satisfied by multi-centered black hole solutions 
of jV = 2 supergravity [24]. This Coulomb branch description is valid provided the centers are 
kept far apart. 

The Higgs branch description, on the other hand, is valid for large values of the chiral multiplet 
scalars ip^, such that the vector multiplets can be integrated out. The moduli space of SUSY vacua 
on the Higgs branch j&h is the set of solutions of the F-term d^W = and D-term equations 

W : £ </>;, T a fe - £ fe T" $ kLa ,, t = a Tr (T«) (4 2) 

Ya>0 Ykt>0 

modulo the action of Yl(U(N(). Equivalently, j$n is the space of semi-stable solutions of d^W = 
modulo Y\(GL(N(,C). BPS states on the Higgs branch correspond to cohomology classes in 
//*(^ H ,Z). The refined index on the Higgs branch is given by the Poincare polynomial (more 
accurately, a polynomial in y, l/y, symmetric under y — > \/y) 

2d 

Q{JZ H \y) = Tr'(-y) 2J i = £ b p {J? H ) (-y)^ . (4.3) 

P =\ 

Here d is the complex dimension of ^h, given by 

d= £ Ytk N f N k - £iV, 2 - / + 1 , (4.4) 

7«>0 t 

where / is the number of independent F-term conditions (/ = for quivers without loop). Since 
./#// is a Kahler manifold, the total cohomology B*[^M^,'£) admits an action of SU(2) known as 
the Lefschetz action, 

J + -h = (0 K /\h , J-=(0 K ^h, J 3 -h= j(n-d)h . (4.5) 

where (Ok is the Kahler form on j$n and n is the form degree of h. The parameter y in (4.3) 
is conjugate to 2/3, and therefore Q{J%H',y) must decompose as a sum of characters of spin j 
representations of SU(2) with j < d/2. 

To demystify these definitions, consider the simplest example of the 2-node quiver with k 
arrows (also known as the generalized Kronecker quiver), with dimension vector (1,1): 

1 k ^ 1 (4.6) 

The vacuum moduli space ^h is the hypersurface £f_j |<^| 2 = c\ in C k , modded out by U(l) 
rotations §t i->- e lQ §i. For c\ < 0, ^u is empty, whereas for c\ > 0, ^u = C^/C x = P^ 1 . The 
complex projective space P k ~ l has bj = 1 for i even, < i < 2(k — 1) and b\ = for i odd. Its 
Poincare polynomial is therefore 

Q{^H\y) = { y> ), y) , (4.7) 

y-l/y 

which is recognized as the character of a spin j = (k— l)/2 representation of the Lefschetz SU(2) 
action. The fact that (4.7) agrees with the Coulomb index (2.23) for two centers is not a coincidence, 
as we shall see shortly. For k = 2 and dimension vector (M,N), one can show that Q(^H',y) is the 
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number of BPS states of charge Mjy +Nj2 in the weak coupling chamber of pure St/ (2) Seiberg- 
Witten theory, again not a coincidence ! 

In order to apply the Coulomb branch formula (3.6) to quiver moduli spaces, we just need 
to identify the Poincare polynomial Q(^H',y) with the total index Q.(y;y), where the charge y is 
identified with the dimension vector {N?} in the positive cone T = N K of the charge lattice Z*\ 
Eq. (3.6) then predicts that the Poincare polynomial Q{^H\y) should be expressible in terms of 
'single centered invariants' Q.s(a), associated to quivers with the same topology but with smaller 
dimension vectors a = {M(}, with M( < Nc for all I = 1, . . .K. While the total indices Q.(y;y) can 
always be traded for the 'single centered invariants' Q.s(a) using (3.6), it is highly non-trivial that 
the £2s(o0's are independent both of the FI parameters q and of the Lefschetz parameter y. In 
the rest of this lecture we shall test this conjecture on a variety of quivers for which Q(^H',y) is 
computable. 

4.2 Quivers without loops 

For quivers without loops, and arbitrary dimension vector, the Poincare polynomial can be 
computed recursively using Harder-Narasimhan nitrations, and counting of points over finite fields 
[40]. This sophisticated mathematics leads to a very explicit formula which can be stated as follows 
[35]: 

1 m r\ 

Qmy,{c&,y) =1 1 -7^vi^^Y\G^{{{Mfy,{c&,y), 



m(y-l/y)^fJr ni8SSKl1 " J ' l ' J "" (4-8) 

E^JLjiK') =#;#(') ||#for i=l,.. .,m 

where the sum runs over all ordered partitions oiN = {N\ , • • • Ng) into parallel vectors M*- ' ' , ■ ■ ■ M^' . 
On the r.h.s, GHiggs({{-^ }; {q }',y) is the 'stack invariant' associated to the quiver with dimension 
vector AfW, given by 

Gm m m};{ce};y) = (- y )-UjWj™*(njfi)-^M 

x (y2_iy-LiM t £ (_l)*-yi^»Ey««K(^.0)WfwyTT(py« )y ]j)-l - (4.9) 

partitions a,i 

Here the sum runs over all ordered partitions of the vector (Mi, ■■■Mi,) into non-zero vectors 
{(N c {, ■ ■ •#£), a = 1, . . . ,s} for s = 1, . . . ,LM, satisfying Nf > 0, ^Nf = M t and 

b K 

££A?C>0 (4.10) 

a=\ i=\ 

for all b between 1 and s — 1. [N,y]l denotes the g-deformed factorial, 

[N,y]\ = [l,y][2,y]...[N,y], [N,y ] = y -^—- . (4.11) 

y L — 1 

When the charge vector {A^} is primitive, Q{{N(};{cf\;y) coincides with GHig g s({^};{Q};>') 
and the corresponding Q ({A^}; {cf\;y) is indeed a symmetric polynomial in v, 1/v with alternating 
sign integer coefficients. If {A^} is not primitive, then the quiver moduli space ^M\\ is singular, and 
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it is more tricky to define its cohomology, nevertheless (4.8) produces a bona-fide Poincare polyno- 
mial in all cases. The Mathematica package supplied with [35] gives a convenient implementation 
of the 'Higgs branch formula' (4.8). 

In order to compare the exact answer (4.8) with the prediction from the Coulomb branch 
formula (3.6), a crucial fact is that the stack invariants (4.9) with arbitrary dimensions {Mg} can all 
be computed in terms of stack invariants of Abelian quivers, a property known as the Abelianization 
(or MPS) formula [15, 17] 

<W{JVi};M;{&};y) = I G Hlggs ({i},{(^ };{(*£/)*' hy)UU-w LL-^ 

{k (t) } i=i i k\'\ \ l \y y ) 

(4.12) 
The quivers appearing on the r.h.s. are obtained by applying the following prescription on each 
node of the original quiver: 

• Each vertex i with dimension Nj is replaced by a collection ig^ of vertices with k = 1 , . . . kg 
for any £ in the partition iV; = Yd^kg 

• Each arrow i — > j (resp. j — > i) with j ^ i in the original quiver induces £ arrows ig^ — > j 
(resp. j — > ig^) in the new quiver for every £,k; 

• The FI parameter at the nodes ig^ is equal to £cu 

Using this formula, one can prove that the Higgs branch formula (4.8) agrees with the Coulomb 
branch formula (3.6) provided the only non-vanishing single-centered invariants £2s(/) w& those 
associated to the basis vectors Jg, with dimension 1 at node £ and dimension at the other nodes. 
These single-centered invariants are moreover moduli-independent and carry zero angular momen- 
tum, namely 

Qs(%;y) = l, ^s(7)=0 if y${y t ,£=l...K}. (4.13) 

Thus, the complete BPS spectrum on the Higgs branch of quiver quantum mechanics associated 
to quiver without loops can be understood from the Coulomb branch picture, where all the states 
are bound states of single-centered black holes carrying unit BPS index. For Abelian quivers, 
this equivalence is expected since the Coulomb branch and the Higgs branch are both regular and 
disconnected, and have overlapping regime of validity [24]. The simplest case of this equivalence 
is the Abelian Kronecker quiver (4.6), with Higgs branch P* _1 , whose de Rham cohomology is 
isomorphic to that the Dirac cohomology of S 2 with k units of flux. More generally, one can show 
that the stack invariant (4.9) for Abelian quivers agrees with the Coulomb index (2.24). The fact 
that this equivalence holds also for non- Abelian quivers is more striking, as the Higgs branch and 
Coulomb branch are no longer disconnected. 

4.3 Quivers with loops 

For general quivers with closed oriented loops and generic superpotential, there is (to my 
knowledge) no general way of computing the cohomology of the quiver moduli space. For Abelian 
quivers however, the situation is much simpler, since the quiver moduli space ^h arises as a 
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complete intersection (corresponding to the F-term constraints) in a product of projective spaces 
^amb (representing the solution of the D-term constraints modulo gauge symmetries). In such 
cases, the Poincare polynomial can be computed by combining the Lefschetz hyperplane theorem, 
which determines all Betti numbers except the middle one, 

b p {J? H ) = l . ' (4.14) 

\bu-p\-Mwab) P>d, 

and the Riemann-Roch theorem to compute %{^h) and hence the middle cohomology bd{^H)- 
Using these techniques, the cohomology of the simplest Abelian quiver with loop 

i (4.15) 

i \ 

c a 

/ X 

1 -*-b — 1 

was studied in [2] and revisited in [25] . Denoting by ^ a ,b,c the quiver moduli space in the chamber 
c\,C2 > 0,C3 < 0, it was found that the generating function of the Poincare polynomials 

Q(x u x 2 ,xy,y) = £ (x 1 ) a (x 2 ) b (x 3 ) c Q(^ aAc ,y) (4.16) 

a>Q,b>Q,c>0 

decomposes into a sum Q = Qc + Qs, where 

xix 2 { 1 -xix 2 +xix 2 xt, (x\ +x 2 +y +y~ l ) } 



Qc{xi,xi,xy,y) 



(1 -x\x 2 )(\ -*i* 3 )(l -x 2 x 3 )(l +xi/y)(l +xiy)(l +x 2 /y){\ +x 2 y) ' 

„2„2„2 

Q&{x\,X2,xy,y) 



(4.17) 



x l x 2 x 3 



(1 — XlX2)(l —X2Xj)(\ — X[X^)[l —X\X 2 —X 2 Xt, —X[Xt, —2x\x 2 x 3 ] 



The first term Qc can be understood as the partition function of the Coulomb indices (2.22) of 3- 
centered black holes carrying charges /i , y 2 , 73 corresponding to the nodes of the quiver, with unit 
degeneracy Q.s{yf,y) = 1. This contribution is both y dependent and moduli dependent (as manifest 
from the lack of invariance under circular permutations of X\,x 2 ,x^), and its Taylor coefficients 
grow polynomially with a,b,c. The second term <2s , on the other hand, is ^-independent, therefore 
associated to Lefschetz singlets in the middle cohomology of ^ a .b,c- Remarkably, it is invariant 
under circular permutations of X\,X2,x$, which shows that the'pure-Higgs' states are robust under 
wall-crossing. They can be shown to contribute only when the triangular inequalities a <b + c,b < 
a + c,c < a + b are satisfied, that is the region where scaling solutions occur. Moreover, their 
contribution to Q(^ a ,b c j) grows exponentially as 2 a+b+c under overall rescaling of a,b,c. The 
total result Q = Qc + Qs agrees precisely with the prediction of the Coulomb branch formula (3.6), 

£(n + Y2 + yny) = gcouiomb(yi,^,^;y) +#({71, n, 73};{i, 1, i};j) +^s(n + 72 + 73) (4.18) 

provided H is chosen according to the minimal modification hypothesis (3.8), and the single- 
centered index Qs(7i + 72 + 73) is equated with the number of Lefschetz singlets. This study was 
extended to all Abelian 'necklace' quivers in [26, 41, 27] with similar results. 

In [27] we have tested the formula on a variety of Abelian quivers with loops, as well as on 
some non-Abelian quivers, and we have found that the Coulomb branch formula (3.6) correctly 
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predicts the cohomology of the quiver moduli space, with constant, ^-independent elementary or 
'pure-Higgs' invariants £2s(/) associated to scaling subquivers. It is worth stressing that unlike 
the case of necklace quivers, the invariants Qs(y) m general do not only contribute to the middle 
cohomology, as they occur as coefficients of non-trivial SU(2) characters in (3.6). 

5. Conclusion and open problems 

In these lectures, we have seen that the simple physical picture of BPS states as bound states 
of elementary constituents carrying their own degrees of freedom and interacting via the super- 
symmetric quantum mechanics of charged particles is remarkably powerful: it provides a simple 
derivation of the wall-crossing formula, and allows to parametrize the BPS index £l(y,t) for all 
values of the parameters or moduli in terms of certain chamber-independent, elementary degenera- 
cies Hs(y). Similarly the protected spin character in JV = 2 gauge theories, or the refined index 
Q.(j;t',y) in jV = 2 string vacua, can be parametrized in terms of refined elementary degeneracies 
Q.s (/;y). While the refined index can get contributions both from short and long multiplets, it seems 
plausible that it should satisfy the same wall-crossing properties as the protected spin character as 
a function of the vector multiplet moduli. On the other hand, it is quite possible that the elementary 
degeneracies £2s(y;y) may jump on certain loci in hypermultiplet moduli space, in particular in 
interpolating from small to large string coupling (at weak coupling, these jumps can only occur in 
co-dimension 2, since the hypermultiplet moduli only affect the holomorphic superpotential W). 
The large string coupling regime is the one where a gravitational picture may apply, and where a 
comparison with the quantum entropy function may be carried out. Unfortunately, the Coulomb 
branch formula (3.6) remains rather impractical, as it requires an enumeration of all possible bound 
states with regular metric. 

In the case of BPS states described at weak coupling by quiver quantum mechanics, the 
Coulomb branch formula is completely well-defined, albeit still conjectural. In all examples that 
we have analyzed, it is satisfying to see that the the elementary degeneracies Cls(j;y) computed 
from the Coulomb branch formula are independent of y, in agreement with the large coupling pre- 
diction that single-centered BPS black holes have zero angular momentum. From a mathematical 
point of view, this is a very unexpected property of quiver moduli spaces, whose understanding may 
help unravel new algebraic structures. From a physics point of view, it would be very interesting to 
clarify the relation between the Higgs and Coulomb branch descriptions and derive the Coulomb 
branch formula (3.6) from quiver quantum mechanics, as it would provide a solvable example of 
open/closed string duality. 

On a longer time-scale, it will be exciting to apply these techniques to actual jV = 2 string 
vacua, such as type II strings compactified on Calabi-Yau threefold, and test agreement between 
macroscopic and microscopic descriptions of BPS black holes, to the same level of precision which 
was achieved for ,jV > 4 string vacua, and possibly address non-BPS black holes as well. 
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